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ABSTRACT 

We study the distribution of angular positions and angular separations of unbound hypervelocity 
stars (HVSs). HVSs are spatially anisotropic at the 3-a level. The spatial anisotropy is significant in 
Galactic longitude, not in latitude, and the inclusion of lower velocity, possibly bound HVSs reduces 
the significance of the anisotropy. We discuss how the observed distribution of HVSs may be linked 
to their origin. In the future, measuring the distribution of HVSs in the southern sky will provide 
additional constraints on the spatial anisotropy and the origin of HVSs. 

Subject headings: Galaxy: bulge — Galaxy: halo — Galaxy: kinematics and dynamics — Galaxy: 
structure — stars: early-type 
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1. INTRODUCTION 



Unbound HVSs were predicted by iHillsl (|1988f ) as the 
natural consequence of the massive black hole (MBH) 
in the Gal actic center. Follo wing the discovery of the 
first HVS (H rown et al]l2005( ). observers have reported 
the discovery of at least 16 unbound HVSs and evi- 
dence for a similar number of bound HVSs ejected by 
the same mechanis m (iHirsch et al.ll2005HEdelmann et al.l 
[2001 iBrown et all l2006aU bl. I2007allbl . \200§) '. Follow- 
up observations of 4 HVSs establish they are main 
seque nce B stars (Fuentcs ct al. 2006\ Bonan os et al.l 
20081: [Przvbilla et al.ll2008bl l(3: lL6pez-Morales fc Bonanod 



200S ) like the S-stars orbiting Sgr A* today ( Ghez et all 



20031 lEisenhauer et al 



ir bitmg 



Martins et al.l 



20081 ). Al- 



though not all unbound stars are necessarily HVSs 
the massive B star HD 271791 is the first example of 
an ur ibound "hyp er-run away" ejected f rom the outer 
disk (jHeber et all 120081 : iPrzvbilia etall l2008ar ) - run- 
away ejection velocities are limited to '^300 k i n s~^ for 
3 Mq stars (Leonard & Duncan' '1988', '1990'; 'Leonard 
1991. 1993; Portcgics Zwart 2000; Davies ct al. 2003 
Guaiandris et all 120051) . Thus the 14 unbound 2.5-4 



M0 stars found in the iBrown et all (|2007bl I2008D tar- 
geted surveys are almost certainly HVSs ejected from 
the Galactic center. 

Remarkably , 8 of the 14 HVSs in the IBrown et al.l 
(|2007U ^Om) targeted surveys are located in just two 
constellations, Leo and Sextans, even though the sur- 
veys cover 1/5*'' of the sky. This spatial anisotropy is 
almost certainly linked to the origin of the HVSs. 

In §2 we show that the observed distribution of HVSs 
on the sky is anisotropic at the 3-<t level. In §3 we dis- 
cuss plausible explanations for the observed anisotropy 
of HVSs. 

2. OBSERVED ANISOTROPY 
2.1. Sample 

We consider the 14 HVSs from the combined surveys of 
IBrown et all ()2007bl l2008t) . Our surveys use the MMT 
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telescope to measure radial velocities for stars with the 
colors of 2.5-4 M© stars. Heliocentric velocities are con- 
verted to Galactocentric velocities assuming that the lo- 
cal rotation speed is 220 km s""' and that the Sun moves 
at (U,V,W)=(10 5.2, 7 .2) km s"' relative to the local 
standard of r est (iDehnen & Binnev 1998). The original 
HVS survey (jBrown et al.ii2007b ) is 100% complete for 
stars with 17 < q'q < 19.5 over 7300 deg^ covered by the 
Sloan Dig ital Sky Survey Data Release 6. The new HVS 
survey (jBrown et al.ll2008[) is 59% complete for stars with 
19.5 < g'g < 20.5 over the same region of sky. 

Figure [T] plots the spatial distribution of stars with 
observed velocities. The combined HVS survey contains 
693 stars and 14 HVSs in the 7300 deg^ S loan region cov- 
ering the north Galactic cap. HVS2 fflirsch et al.ll2005h 
is also located in this region (see Fig. |T]), however it falls 
outside our color / magnitude criteria. Thus we exclude 
HVS2 from this analysis. 

2.2. Significant Anisotropy 

Figure [2] plots the cumulative Galactic longitude and 
latitude distributions of the HVSs and the other survey 
stars. Kolmogorov-Smirnov (K-S) tests find 0.007 and 
0.11 likelihoods that the HVSs are drawn from the same 
longitude and latitude distributions, respectively, as the 
survey stars. Thus the distribution of HVS longitudes 
appears anisotropic at the 3-cr level. 

As a second test, we explore the anisotropy in terms of 
the distribution of angular separations, 9, of the HVSs 
compared to the survey stars. Because the new HVS sur- 
vey is not yet complete, we calculate 0's for all unique 
pairs of stars in the new and original surveys separately. 
The original survey includes HVS4-HVS10; a K-S test 
finds a 0.031 likelihood that those HVSs are drawn from 
the same distribution of 9 as the original survey stars. 
The new survey includes HVSl and HVS11-HVS16; a K- 
S test finds a 7 x 10~^ likelihood that those HVSs are 
drawn from the same distribution of 9 as the new survey 
stars. Figure [3] plots the 9^s of both surveys concatenated 
together. The likelihood of the combined set of HVSs is 
7 X 10^^; thus the distribution of HVS angular separa- 



2 



Brown et al. 



0" 




J = 180" 



Fig. 1. — Polar projection, in Galactic coordinates, showing 
the 14 unbound HVSs (stars) and the 693 other stars (diamonds) 
in our HVS survey (iBrown et al.ii2008i) covering the north Galactic 
cap. HVS2, while not part of the survey, is also marked (plus sign). 
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Angular Separation 9 [deg] 

Fig. 3. — Upper panel: cumulative distribution of angular sep- 
arations, 9. Lower panel: two-point angular correlation function 
w(6) of the HVSs with respect to sets randomly drawn from the 
survey region. 
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Fig. 2. — Cumulative distributions of Galactic I and b of the 
14 HVSs (da shed lines) and the 693 other stars (solid lines) in our 
HVS survey IIBrown et al.ll2008lV 



tions differs from the distribution of survey star angular 
separations at tlie 5-a level. 

As a third test, we measure the clustering of HVSs us- 
ing the two-point angular correlation function iu{9). We 
use a Monte Carlo estimator (Landy & Szalay 1993j) and 
compare the observed HVSs against 10^ sets randomly 
drawn from the survey region. The lower panel of Fig- 
ure [3] plots the resulting w{9) in 15° bins. Errorbars are 
determined by Poisson statistics. HVSs are clustered at 
small angular separations 9 < 45° and missing at large 
angular separations 9 > 60° with ^3.5-(t significance. 

2.3. Velocity Dependence 

We now consider the spatial aniso tropy of lower veloc- 
ity stars that may also be HVSs. iBrown et al.1 ()2008D 
identify 4 "possible HVSs," stars that are bound in the 
Kenvon et al. (2008') potential model but unbound in the 
Xue et al. (2008) potential model. Adding the 4 possible 
HVSs to the above analysis reduces the significance of 
the anisotropy to the 2-a level. There are also 8 possi- 
bly "bound HVSs," stars with Vrf > +275 km s~^ that 
are s ignificant o utliers from the overall velocity distribu- 
tion (jBrown et al. 2008). Adding the bound HVSs to the 
above analysis yields an insignificant anisotropy. Thus 
lower velocity stars have a more isotropic distr ibution, a 
trend noted previously in IBrown et al.l (|2007a[ ). 



2.4. HVS Pairs 

There are 3 pairs of unbound HVSs with angular sep- 
arations less than 3.5° (see Fig. H]): HVS7 & HVS15 
near (;,6) = (265°, 55°), HVS12 & HVS13 and HVS12 & 
HVS14 near (/, b)= (245°, 52°). Any physical association 
between the individual HVSs, however, appears unlikely. 
HVS7 & HVS15 are separated by 2.5° but have velocities 
and distances that imply a ~70 Myr differ ence in travel 
time from the Galactic Center (see Fig. 3 of lBrown et al.l 
[20081 ). HVS12 has a 429 km s^^ minimum rest-frame ve- 
locity very similar to that of HVS13 and HVS14, but it 
has half the distance of the other two HVSs. Thus none 
of the HVS pairs share a common ejection event. 

3. ORIGIN OF THE SPATIAL ANISOTROPY 

We observe that the spatial anisotropy of unbound 
HVSs is statistically robust, that lower velocity HVSs are 
systematically more isotropic, and that apparent close 
pairs of HVSs are physically unrelated. Possible expla- 
nations for the observations include: 
Selection Effect. Previously, we argued that the HVS 
anisotropy may be a selection effect of our magnitude- 
limited survey and the Sun 's off-center location in the 
Galaxy (jBrown et al.ll2007al ). However, this selection ef- 
fect can account only for an extra ~10% HVSs in the 
anti-center hemisphere, not all of the HVSs in the anti- 
center hemisphere. IVIoreover, the observed HVSs cluster 
around I ^ 240°, not I = 180°. 

Runaways. Runaway stars like HD 271791 may con- 
taminate the population of HVSs. However, we expect 
runaways ejected from the disk to have an isotropic dis- 
tribution in G alactic longitude, as demonstrated by the 
llVIartinI ()2006( ) Hipparcos-selected sample of runaway B 
stars. Moreover, because r unaways are systematically 
ejected at low velocities (e.g.. lPortegies Zw art 2000), the 
fastest runaways are those ejected in the direction of 
Galactic rotation and thus preferentially found at low 
Galactic latitudes. Thus the expected distribution of 
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runaway longitudes and latitudes are contrary to the ob- 
served distribution of HVSs. 

Large Scale Structure. The distribution of Local Group 
dwarf gal axies is anisotropi c, possibly due to a tidal ori- 
gin (e.g., iMetz et al] I2008D . A tidal debris origin ap- 
pears supported by the clumping of HVS travel times 
around 100-200 Myr, however the travel times are sim- 
ply a product of the HVS's '^500 km s~^ velocities and 
our magnitude-limited survey depth of 50-100 kpc. HVS 
travel times are in fact problematic for a tidal debris ori- 
gin because the times are a significant fraction of the 
stars' main sequence lifetimes, and multiple (gas-rich) 
tidal disruption events would be required to explain the 
full 2 X 10* yr span of HVS travel times. No dwarf galaxy 
in the Local Group travels with radial velocities compara- 
ble to the unbound H VSs; known dwarf galaxy remnants 
hke the Sgr stream (jibata et al.lll994f ) are bound. We 
thus consider tidal debris an unlikel y explanat ion for the 
observed set of HVSs (however, see lAbadi et a l. 2008). 
Binary Black Hole. While an equal -mass binary MBH 
is rul ed out in the Galactic Center (jReid &: Brunthaleil 
l2004f ). theorists speculate that the massive star clus- 
ters in the Galactic Center form intermediate mass black 
holes (IMBHs) in their cores. If such IMBHs exist, dy- 
namical friction causes them to in-spiral into the cen- 
tral MBH, preferentially ejecting HVSs from their or- 
bital planes. Thus the expected signat ure of a IMBH in- 
spira l is a r ing of HVSs around the sky ("Gualandri s et all 
120051: iLevin 2006; Sesana et al. 2006). Baumgard t" et all 
1(2001 argue, however, that stellar interactions perturb 
the orbital plane of an in-spiraling IMBH; the result- 
ing HVS distribution in this scenario may in fact be 
isotropic. Moreover, a single IMBH in-spiral event hap- 
pens on timescales 10-100 x shorter than the observed 
span of HVS travel times; multiple IMBH in-spiral events 
are required to explain the observed HVSs. 
Galactic Center Structure. The Galactic center con- 
tains many we l l-defin ed structures. As illustrated in 
iPaumard et al.l (|2006f ). the molecular gas circum- nuclear 
disk and the ionized northern arm are roughly aligned 
with the plane of the Milky Way. The gaseous minispiral 



is perpendicular to the plane of the Milky Way. Notably, 
the stellar disk 0.1 pc from the MB H is roughly p erpen- 
dicular to the gaseous components (iLu et al.ll2008 ) . The 



stellar disk contains massive stars ( Tanner et al 
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IPaumard et al.l[2006D. possibly formed in-situ from a gas 
accre tion disk (jGenzel et al.ll2003l : iLevin fc Beloborodovl 
l2003f ). Dynamical interactions between a pair of stellar 
disks may scatter stars in toward s the MBH, explaining 
both the S-stars and the HVSs (|Lockmann et al.l 120081 : 
iPerets et al.ll200 8^. Clearly, the Galactic center contains 
non-isotropic distributions of stars and gas which may 
provide a natural source for the observed anisotropy of 
HVSs ejected from the Galactic center. However, it is un- 
clear if the observed structures can persist long enough 
to explain the anisotropic distribution of HVSs. 

4. CONCLUSION 

Unbound HVSs are spatially anisotropic at the 3-cr 
level. The anisotropy is most significant in Galactic lon- 
gitude, and not in latitude. Lower velocity HVSs are 
systematically more isotropic, and apparent close pairs 
of HVSs are physically unrelated. 

The observed d istribution o f HVSs is linked to the ori- 
gin of the HVSs. lAbadi et all (|2008l ) propose a tidal de- 
bris explanation, although this appears difficult to rec- 
oncile with all the observations. We investigate other 
physical models for the anisotropy in a separate paper. 
In the future, measuring the distribution of bound and 
unbound HVSs over the southern sky will allow us to 
better constrain the anisotropy and the origin of HVSs. 



This work is based on observations obtained at the 
MMT Observatory, a joint facility of the Smithsonian 
Institution and the University of Arizona. This research 
makes use of NASA's Astrophysics Data System Biblio- 
graphic Services. We thank the anonymous referee and 
Gleg Gnedin for helpful comments. This work was sup- 
ported by the Smithsonian Institution. 

Facilities: MMT (Blue Channel Spectrograph) 



REFERENCES 



Abadi, M. C, Navarro, J. F., & Steinmetz, M. 2008, ApJ, 
submitted 

Baumgardt, H., Gualandris, A., & Portegies Zwart, S. 2006, 

MNRAS, 372, 174 
Bonanos, A. Z., Lopez-Morales, M., Hunter, L, & Ryans, R. S. L 

2008, ApJ, 675, L77 
Brown, W. R., Geller, M. J., & Kenyon, S. J. 2008, ApJ, accepted 
Brown, W. R., Geller, M. J., Kenyon, S. J., & Kurtz, M. J. 2005, 

ApJ, 622, L33 
— . 2006a, ApJ, 640, L35 
— . 2006b, ApJ, 647, 303 

Brown, W. R., Geller, M. J., Kenyon, S. J., Kurtz, M. J., & 

Bromley, B. C. 2007a, ApJ, 660, 311 
— . 2007b, ApJ, 671, 1708 

Davies, M. B., King, A., & Ritter, H. 2002, MNRAS, 333, 463 

Dehnen, W. & Binney, J. J. 1998, MNRAS, 298, 387 

Edelmann, H., Napiwotzki, R., Heber, U., Christlieb, N., & 

Reimers, D. 2005, ApJ, 634, L181 
Eisenhauer, F. et al. 2005, ApJ, 628, 246 

Fuentes, C. L, Stanek, K. Z., Gaudi, B. S., McLeod, B. A., 
Bogdanov, S., Hartman, J. D., Hickox, R. C, & Holman, M. J. 
2006, ApJ, 636, L37 

Genzel, R. et al. 2003, ApJ, 594, 812 

Chez, A. M. et al. 2003, ApJ, 586, L127 



Gualandris, A., Portegies Zwart, S. P., & Sipior, M. S. 2005, 

MNRAS, 363, 223 
Heber, U., Edelmann, H., Napiwotzki, R., Altmann, M., & Scholz, 

R.-D. 2008, A&A, 483, L21 
Hills, J. G. 1988, Nature, 331, 687 

Hirsch, H. A., Heber, U., O'Toole, S. J., & Bresolin, F. 2005, A&A, 
444, L61 

Ibata, R. A., Gilmore, G., & Irwin, M. J. 1994, Nature, 370, 194 
Kenyon, S. J., Bromley, B. C, Geller, M. J., & Brown, W. R. 2008, 

ApJ, 680, 312 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 
Leonard, P. J. T. 1991, AJ, 101, 562 

Leonard, P. J. T. 1993, in ASP Conf. Ser. 45, Luminous High- 
Latitude Stars, ed. D. Sasselov, 360 
Leonard, P. J. T. & Duncan, M. J. 1988, AJ, 96, 222 
— . 1990, AJ, 99, 608 
Levin, Y. 2006, ApJ, 653, 1203 

Levin, Y. & Beloborodov, A. M. 2003, ApJ, 590, L33 

Lockmann, U., Baumgardt, H., & Kroupa, P. 2008, ApJ, 683, L151 

Lopez-Morales, M. & Bonanos, A. Z. 2008, ApJ, 685, L47 

Lu, J. R., Chez, A. M., Hornstein, S. D., Morris, M. R., Becklin, 

E. E., & Matthews, K. 2008, ApJ, accepted 
Martin, J. C. 2006, AJ, 131, 3047 

Martins, F., Gillessen, S., Eisenhauer, F., Genzel, R., Ott, T., & 
Trippe, S. 2008, ApJ, 672, L119 



4 



Brown et al. 



Metz, M., Kroupa, P., & Libeskind, N. I. 2008, ApJ, 680, 287 
Paumaxd, T. et al. 2006, ApJ, 643, 1011 

Perets, H. B., Gualandris, A., Merritt, D., & Alexander, T. 2008, 

astro-ph/0807.2340 
Portegies Zwart, S. F. 2000, ApJ, 544, 437 

Przybilla, N., Nieva, M. F., Heber, U., & Butler, K. 2008a, ApJ, 
684, L103 

Przybilla, N., Nieva, M. F., Heber, U., Firnstein, M., Butler, K., 
Napiwotzki, R., & Edelmann, H. 2008b, A&A, 480, L37 



Przybilla, N., Nieva, M. P., Tillich, A., Heber, U., Butler, K., 

Brown, W. R. 2008c, A&A, 488, L51 
Reid, M. J. & Brunthaler, A. 2004, ApJ, 616, 872 
Sesana, A., Haardt, F., & Madau, P. 2006, ApJ, 651, 392 
Tanner, A. et al. 2006, ApJ, 641, 891 
Xue, X. et al. 2008, ApJ, 684, 1143 



